Quantum versus classical behavior in the boundary susceptibility of the ferromagnetic 

Heisenberg chain 

J. SirkeiQ 

Department of Physics and Astronomy, University of British Columbia, Vancouver, B.C., Canada V6T 1Z1 



M. Bortz 

Department of Theoretical Physics, Research School of Physics and Engineering, 
Australian National University, Canberra ACT 0200, Australia 
(Dated: February 2, 2008) 

We calculate the temperature dependence of the boundary susceptibility \b for the quantum 
ferromagnetic Heisenberg chain by a modified spin- wave theory (MSWT). We find that \b diverges 
at low temperatures ~ — T~ 3 and therefore more rapidly and with opposite sign than the bulk 
susceptibility xtmik ~ T~ 2 . Our result for \b is identical in leading order with the result for the 
classical system. In next leading orders, however, quantum corrections to the classical result exist 
which are important to obtain a good description over a wide temperature range. For the S — 1/2 
case, we show that our full result from MSWT is in excellent agreement with numerical data obtained 
by the density-matrix renormalization group applied to transfer matrices. Finally, we discuss the 
quantum to classical crossover as well as consequences of our results for experiment in some detail. 
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I. INTRODUCTION 

Although the classical and the quantum version of the 
ferromagnetic Heisenberg model have the same ground 
state, it is not obvious in how far the low temperature 
properties of these systems are also similar. The Hamil- 
tonian for the quantum ferromagnetic chain with open 
boundary conditions (OBC), spin S and N sites in a 
magnetic field h along the z-direction is given by 
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H = - J S n S n+ i - g ^2 hS. 
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where J > is the coupling constant and g the g-factor. 
The classical version of this model can be obtained by in- 
troducing new unit vector operators s n = S n /S^ These 
new operators commute in the limit S — > oo, leading to 
the classical Hamiltonian 
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To allow for a comparison between the quantum and the 
classical Hamiltonian for different spin values S we have 
set J c = JS 2 and g c = gS. 

The finite temperature properties of the classical model 
have been calculated several decades ago for OBC— as 
well as periodic boundary conditions (PBC) 2 . Whereas 
the first correction to the total free energy in the ther- 
modynamic limit for PBC is 0(1/N), the free energy for 
OBC contains a term 0(1). This boundary or surface 
free energy then yields 0(l)-contributions to all other 
thermodynamic quantities as for example the suscepti- 
bility. From Fisher's results 1 - one finds that the classical 



bulk susceptibility behaves as 
Xbuik = N 9. 
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whereas the classical boundary susceptibility (0(1) con- 
tribution) is given by 
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Xb therefore diverges more rapidly and with opposite 
sign than x£ ulk - 

A very different behavior for bulk and boundary 
susceptibility has also recently been observed for the 
quantum antiferromagnetic S = 1/2 XXZ-chain with 
anisotropy < A < l^dik For this system it is known 
that the bulk susceptibility is finite for T — > with the 
T = value of Xbuik depending on the anisotropy A. 
Xb(T), on the other hand, is finite only for < A < 1/2 
whereas it diverges for 1/2 < A < 1 when T — > 0. By 
a combination of different techniques like bosonization, 
conformal field theory, Bethe ansatz as well as numerical 
results, a complete picture of the low-temperature prop- 
erties of xb has been obtained. 3 ! 4 ! 5 These results are not 
only of theoretical interest but might also be relevant 
for realizations of quasi one-dimensional antiferromag- 
nets as for example SrCuOa, in particular, when such 
compounds are doped with a moderate amount of non- 
magnetic impurities. In such a case the spin chain will be 
partitioned into finite chains with essentially free bound- 
aries and knowledge of the boundary contributions will 
be essential to understand experiments on such systems. 

For the quantum S — 1/2 ferromagnetic chain the 
standard bosonization approach and conformal field the- 
ory are not applicable because the dispersion relation is 
quadratic instead of linear. The model is, however, still 
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integrable and thermodynamic properties can in prin- 
ciple be calculated either by the thermodynamic Bethe 
ansatz 6 (TBA) or by the Bethe ansatz applied to quan- 
tum transfer matrices i The bulk susceptibility has in- 
deed been obtained by an analysis of the TBA equations. 
In the first analysis of this kind by Schlottmann& it has 
been proposed that xtmik ~ J/T 2 \n(J/T). A later nu- 
merical analysis of the TBA&^i^, however, showed that 
Xbuik ~ J/6T 2 at low temperatures as in the classical 
case but with corrections to this leading term which are 
different from Eq. ©. In addition it has been found 
that the leading term as well as the quantum correc- 
tions can be obtained by a modified spin-wave theory 
(MSWT)iS*ii It has later been shown that the classical 
and quantum ferromagnetic chains obey the same scal- 
ing laws at low temperatures. 12 Furthermore, the critical 
theory controlling the low-energy behavior of both chains 
has been identified, which explains in more general terms 
why Xbuik an d Xbuik are identical at low temperatures^ 

It is still unclear how the TBA has to be modified to 
allow also for the calculation of boundary contributions. 
Some of th e difficulties one encounters are discussed in 
Refs. I r ll ll Within the QTM approach an explicit for- 
mula for the boundary free energy has been derived very 
recently»ii The explicit evaluation of this formula, how- 
ever, is still a formidable task because it involves expec- 
tation values of an operator in the dominant eigenstate 
of the QTM which are notoriously difficult to calculate. 

For these reasons we will follow here a different route 
and will use in section [H] Takahashi's MSWT, which 
has been so successful for the bulk, to calculate the 
boundary susceptibility. In section UTTl we then compare 
our result with numerical data obtained by the density- 
matrix renormalization group applied to transfer matri- 
ces (TMRG). In the last section we discuss the quantum 
to classical crossover observed and comment on the rele- 
vance of our results for experiment. 



II. MODIFIED SPIN- WAVE THEORY 

With the help of the Holstein-Primakoff transforma- 
tion 

5+ = V2S\J 1 - aia n /2S a„ , — S - a) n a n (5) 

the model JIJ can be represented exactly in terms of 
bosons a n . Linear spin- wave theory is obtained if one 
replaces the second square root in Eq. JSJ by 1. Correc- 
tions to this simple approximation can be calculated in 
principle in a systematic way by expanding the square 
root in powers of 1/S. In any of these approximations 
it is important to notice that the bosons have to obey 
a hard-core constraint restricting the maximum number 
of bosons per site to 2S. In higher dimension it is often 
acceptable to ignore this constraint completely. In one 
dimension, however, this constraint is crucial but hard 
to incorporate locally. Because the SU (2) symmetry in 



a system with h — can only be broken at T = we 
might instead try to introduce a constraint which fixes 
the number of bosons to be S on the average so that 

n 

at any finite temperature. This approach has been used 
successfully by Takahashi 6 - to calculate the free energy 
and the susceptibility for a chain with PBC. We will use 
the same approach here for a system with OBC to obtain 
the boundary susceptibility. 

Let us first rederive Takahashi's result for PBC in a 
slightly different way. Expanding up to quartic order 
in the boson operators in Eq. (JjJJ and using a one-loop 
approximation for the quartic terms, we obtain for the 
Hamiltonian JQ) at zero magnetic field 

H = JS' t{k)a\a k + v a\a k 

k k 

e{k) = 2(1- cos k) (7) 

where 

k 

and the average number of bosons n k is given by 

n k := (o^Ofc) = [exp(JS'e(fc)/T + v) - l]" 1 . (9) 

At temperatures T/J < 1 the number of bosons in high 
momentum states is small. The bosons in low momentum 
states, on the other hand, will yield only a small contri- 
bution to the sum in Eq. JHJ so that we will set S' = S in 
the following. According to Eq. ©, the potential v then 
has to be determined in such a way that 

S=^$>. (10) 

fc 

Differentiating the partition function for the Hamiltonian 
(JTJ twice with respect to h one finds that the suscepti- 
bility is given by X = .fV^E^^^). However, the 
spin- wave expansion we are using here breaks the SU (2) 
symmetry so that we will calculate the susceptibility in- 
stead by 

2 f N 1 

X = | EE < 5 „S m ) + NS(S + 1) \ . (11) 

I n— 1 m^n J 

In this way the consequences of SU (2) symmetry break- 
ing are less severe due to the averaging over all three 
directions. Using the constraint ifTT)|) one finds 6 
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The momenta for a chain of length N with PBC are given 
by k = 2nl/N where I = 0, 1, • • • , N - 1. For T/J < 1 
the most important contributions to the sum in l|12|) come 
from k « 0, 2tt and we can evaluate these contributions 
by using a saddle point integration 



cos[/c(r„ - r m )] 
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In the last line we have introduced the abbreviation 
t = T / J S. To understand why the saddle point approxi- 
mation for the integrand in (|13J) is indeed sufficient here 
we make the following observation: l|13fl can be evaluated 
alternatively by closing the integration contour in the up- 
per or lower half of the complex A:-plane, depending on 
the sign of r n — r m . The residues closest to the real axis 
then yield <|14|l . Next-leading residues give contributions 
0(VTexp[-VT}). These would result in terms 0(1/VT) 
in the susceptibility, which are neglected in the ongoing. 

From the constraint IjlUfl one can easily determine the 
potential v as a series in y/i. The result is& 

^=!+<(!) J +< 2 (!) 3 +°^ < i5 > 

where q = £(l/2)/y / 7r. When we rewrite the correlation 
function (|14fl in terms of the normalized spin operators s n 
and the coupling constant J c as given in the introduction 
and use only the leading term from 115|) we find 



(s n s m ) = exp(-|r n - r m \T/J c ) 



(16) 



for all values of 5*. In particular, the correlation length 
at T/J <C 1 is always given by £ = J c /T. Further- 
more, Eq. (|16[) also agrees with the result for the classical 
modelti Note, however, that this is no longer the case if 
one takes the next-leading terms in l|15|l into account. 

To calculate the susceptibility we have to evaluate the 
sum in Eq. (|llfl . For PBC each distance \r n — r m \ = 
1, • • ■ , N/2 appears 2N times. The susceptibility in the 
thermodynamic limit can therefore be obtained by 

2 f N/2 



w^oo 3T 1 2v 
Ng 2 
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The first term agrees exactly with the result obtained by 
Takahashi 6 , however, we find here in addition the sec- 
ond and third term, which are absent in Takahashi's re- 
sult. Note, that these terms exactly cancel each other 



for S — > oo. The differences between our and Taka- 
hashi's result can be explained as follows: Whereas in 
Ref. |6j the sum in Eq. (|11|) is carried out without ap- 
proximating the correlation function l|12ll we have taken 
here only the long-distance asymptotics of (S n S m ) into 
account as obtained by the saddle point approximation 
in H14J) . Interestingly, the terms in the susceptibility up 
to 0(1/T) remain unaffected, i.e., these terms are not 
influenced by the behavior of the correlation function 
at short-distances. In fact, we might trust our spine- 
wave approximation only in the long-wavelength limit 
where the spin-wave interaction is small. In our one- 
loop approximation this becomes clear when considering 
Eqs. |0BJ. When all momenta involved are small, the 
sum in Eq. (jHJ is also small and S' rs S. In this limit the 
Hamiltonian becomes equivalent to the one for ideal 
non-interacting spin waves. 

For these reasons we cannot expect that the MSWT 
gives reasonable results if we try to calculate local quan- 
tities for OBC near the boundary. We observed that 
neither a local constraint (S*) = nor the correlation 
function (S n S m ) can be calculated without inconsisten- 
cies. For example, if we calculate the correlation function 
for OBC explicitly we find a constant term which van- 
ishes only if we set v = t/AS 2 exactly. However, the 
condition l|lUfl still requires corrections to v as given in 

03- 

Far enough away from the boundaries, on the other 
hand, the correlation function will still behave as in 
Eq. (13) . When we perform the sum in l|ll(l using again 
this long-distance asymptotics for (S n S m ) but in a way 
appropriate for OBC we will already obtain a 0(1) cor- 
rection to the susceptibility without taken the modifica- 
tions to the correlation function near the boundary into 
account. We conjecture that for low temperatures, this 
term yields xb- The physical picture behind this proce- 
dure is as follows: We can combine two open chains each 
of length M — 1 to one periodic chain of length N = 2M, 
where the two additional sites do not couple with their 
neighbors. We then carry out the sum in Eq. I|ll|) only 
over one half of the periodic chain, thereby discarding 
correlations between this subsystem and the rest. Doing 
so we ignore local differences between PBC and OBC. 

What makes us confident that this is indeed sufficient 
to obtain the leading terms in a low-temperature expan- 
sion for the boundary susceptibility xb is that the leading 
term ~ — 1/T 3 is universal in the sense that it does not 
depend on S. Especially, it is the leading term of Xb 
for both 5 = 1/2 and S = ooi 12 i 13 The classical result 
for the correlation function (|16f) has been first obtained 
by Fisher 1 for an open chain. I.e., in the classical limit 
the exponential decay of the correlation function does de- 
pend only on |r„ — r m \ and not on r n , r m alone, although 
translational invariance is broken! 

We therefore conjecture that the leading terms in a 
low-temperature expansion of the susceptibility for a 
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quantum chain with OBC are given by 



Xobc = -ztp < 
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In particular, the boundary susceptibility is given by 

,2 



XB 



9 



24J5 
, 25 3 



-J-'v- 2 



3Jt 3 



2\/t 2 3* 



(19) 



Note that the leading term is identical to the leading 
term in the classical result (@J when J, g are replaced by 
J C ,5 C . This confirms our expectations. To test if the pro- 
cedure proposed here gives indeed the right corrections 
to the classical result we will check formula ff 1 01) against 
numerical data for the 5 = 1/2 quantum model in the 
following section. 



III. NUMERICAL RESULTS 

In a system with OBC the one-point correlation func- 
tion (5 z (r)) is no longer a constant because translational 
invariance is broken. We define 



C(r) = <S*(r)> 



OBC 



(20) 



where m is the magnetization per site in the system with 
PBC and r is the distance from the boundary. The 
local boundary susceptibility is then given by xb{t) = 
dC(r)/dh\h=o and the total boundary susceptibility xb 
can be obtained by 
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This means that we can calculate xs by considering only 
a local quantity which is particularly useful in numerical 
calculations where it is difficult to obtain the 0(1) contri- 
bution directly. Particularly suited for this purpose is the 
density-matrix renormalization group applied to transfer 
matrices (TMRG) because the thermodynamic limit is 
performed exactly. The idea of the TMRG is to express 
the partition function Z of a one-dimensional quantum 
model by that of an equivalent two-dimensional classi- 
cal model which can be derived by the Trotter-Suzuki 
formulaii£*i& For the classical model a suitable transfer 
matrix T can be defined which allows for the calculation 
of all thermodynamic quantities in the thermodynamic 
limit by considering solely the largest eigenvalue of this 



transfer matrix. Details of the algorithm can be found 
in Refs. Il7lfl3 ll9 20. The method has been extended to 
impurity problems in Ref. |2l[ In particular, the local 
magnetization at a distance r from the boundary of a 
system with N sites can be obtained by 



(22) 



where \^\) ((^£1) are the right (left) eigenstates of the 
transfer matrix T, T is a modified transfer matrix con- 
taining the broken bond and T(5 Z ) is the transfer matrix 
with the operator S z included. Because the spectrum of 
T has a gap between the leading eigenvalue A and the 
next-leading eigenvalues, Eq. H22[> reduces in the thermo- 
dynamic limit to 
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Therefore only the leading eigenvalue and the corre- 
sponding eigenvectors have to be known to calculate the 
local magnetization in the thermodynamic limit. Far 
away from the boundary (S z (r)) becomes a constant, the 
bulk magnetization 



lim lim (S z (r)) 

r— >oo N—^oo 



= lim 
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To obtain numerically the susceptibility profile xb{t) we 
calculate C(r) for small fields h ~ 10~ 4 , 10 -5 by using 
Eqs. (|23I24|I and then taking the numerical derivative. 

Here we want to study the quantum model Q with 
S = 1/2, J = 1 and g = 2. First, we want to test our nu- 
merical results by calculating the bulk susceptibility and 
comparing with Eq. I|17l) which agrees with the TBA. 6 
The result is shown in Fig. ^ and the agreement at low 
temperatures is excellent. Note also, that although the 
leading terms in the low-temperature expansion for the 
classical and the quantum model are identical, extremely 
low temperatures are necessary to see the classical scaling 
for the S =1/2 quantum model. 

The boundary susceptibility is shown in Fig. 0m com- 
parison to the classical result as well as to formula H19|) 
conjectured for the quantum case. The excellent agree- 
ment confirms our conjecture for the 5=1/2 case. As 
()19J) also agrees with the classical result in the limit 
S — > oo we expect that our result is valid for all S. 

Finally, we show in Fig. [3] susceptibility profiles xb{t) 
for different temperatures. As the total boundary sus- 
ceptibility is given by Eq. H21(l which, on the other hand, 
should be equal to i|19|) we can even determine an analytic 
formula for xb (?") and find 
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FIG. 1: Xbuik as a function of temperature. The circles denote 
the numerical data obtained by TMRG, the dashed line is 
the classical result from Eq. and the solid line Takahashi's 
result G3 obtained by MSWT. 




T/J 



FIG. 2: Boundary susceptibility \b as a function of tem- 
perature. The circles denote the numerical data obtained by 
TMRG, the dashed line is the classical result from Eq. (gj 
and the solid line our result H19H from MSWT. 



This formula is in excellent agreement with our numerical 
data (see dashed lines in Fig. • The deviations at large 
distances r where Xb{ t ) is small are due to numerical 
errors. For the fields h ~ used here, xb{t) ~ 10~ 2 
corresponds to a local magnetization C(r) ~ which 
becomes comparable with the accuracy of the calculation. 
Note also, that according to Eq. 1|25[) the one-point cor- 
relation function (S z (r)) will decay for small magnetic 
fields with exactly the same correlation length as the 
bulk two-point correlation function (S z (r)S z (0)) . This 
connection between one and two-point correlation func- 
tion has also been found for the antiferromagnetic XX Z 
chain£ 




r 



FIG. 3: Susceptibility profile Xs(r) at a distance r from the 
boundary for different temperatures T = 0.013, ■ ■ ■ , 0.8. The 
dots represent the numerical, the dashed lines the theoretical 
result according to Eq. I|25|l . 



IV. CONCLUSIONS 

We want to emphasize that the boundary susceptibility 
is not a finite size quantity. It is defined as the difference 
in susceptibilities between a periodic chain and a chain 
with OBC in the thermodynamic limit. In fact, when we 
calculated Xpbc and Xobc in section |n] we have ignored 
terms ~ exp(— 2Ny/tv). For a finite chain with OBC this 
is a valid approximation if T/J > 1/4 AT and our results 
can be directly applied if this condition is fulfilled. At 
temperatures T/J ~ l/N, where finite size corrections 
are sufficiently small to be ignored, we find a ~ 25% re- 
duction of the total susceptibility in the open compared 
to the periodic system. This effect should therefore be 
relevant in susceptibility measurements on systems with 
non-magnetic impurities when the temperature T/J be- 
comes comparable to the concentration of impurities (in- 
verse average chain length). 

In this context we want to mention that the low- 
T behavior of x PB c following from Eqs. I|15I17|) has 
been observed experimentally! 2 ^ 2 ^ Furthermore, con- 
trolled doping of quasi-one dimensional ferromagnets 
with both magnetic^ and non-magnetio2^ defects is pos- 
sible. Most interestingly, susceptibility measurements 
of diluted two-dimensional ferromagnets have revealed a 
one-dimensional behavior at the percolation threshold, 26 
and an unexplained lowering of the susceptibility under 
the percolation threshold at low temperatures. It would 
be certainly interesting to try to understand these exper- 
iments in more detail in the light of the results presented 
here. 

Finally, we want to address the question at which 
temperature scale the crossover from quantum to clas- 
sical behavior occurs. Clearly, the system behaves classi- 
cally at length scales much smaller than the correlation 
length £ = 1/2^/tv w J c /T where all spins are practi- 
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cally aligned. The length scale fo r fluctu ations is set by 
the spin-wave wavelength A ~ y/ J C /TS. So we expect 
classical behavior when A <C £ which is true for all S at 
sufficiently low temperatures. As expected, A becomes 
smaller with increasing S whereas the correlation length 
£ does not change. Therefore the crossover temperature 
will increase with the spin quantum number S. 

In summary, we have used a modified spin-wave theory 
- where a chemical potential guarantees zero magnetiza- 
tion at zero magnetic field for any finite temperature - 
to calculate the boundary susceptibility xb for the open 
spin-5 quantum ferromagnetic chain. We found that xb 
can be expanded in powers of vT and that the lead- 
ing term is given by xb ~ — 1 /T 3 in agreement with the 
classical result. The quantum corrections to this classical 
result are, however, important to obtain a good descrip- 
tion over a wide temperature range. We have verified our 
formula for the 5=1/2 case by comparing with numeri- 



cal data obtained by the density-matrix renormalization 
group applied to transfer matrices and have found ex- 
cellent agreement. We have been even able to derive 
an analytic formula for the local boundary susceptibil- 
ity Xb{t) which we also checked numerically. Most im- 
portant, we have shown that xb at low temperatures is 
"universal" , in the sense that it is completely determined 
by the long-distance asymptotics of the two-point corre- 
lation function (S n S m ). 
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